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RESISTIVE  INTERCHANGE  MODES  IN 
REVERSED  FIELD  PINCHES 


I.  Introduction 

It  has  long  been  recognized  that  reversed  field  pinches,  having  unfavorable  average  curvature 
should  be  unstable  to  resistive  interchange  modes  for  arbitrarily  small  0  and  for  arbitrarily  large  m  and 
n  (poloidal  and  toroidal  mode  numbers)  unless  finite  Larmor  radius  effects  are  taken  into  account.1  2 
The  theory  of  gravity  driven  modes  in  an  incompressible  fluid  in  a  sheared  field  indicates 

fBz 

y  —  i) 1/3 1/? '/VI2'0  where  y  is  the  growth  rate,  tj  the  resistivity,  p'  the  pressure  gradient  and  q  -  ■■  ■  . 

J\Dg 

However  recent  numerical  results  pertaining  to  the  linear  and  nonlinear  behavior  of  these  modes  by 
Schnack  et  al.3  indicate  that  at  low  0,  the  plasma  is  in  fact  much  more  stable  than  this  result  would 
indicate.  We  find  that  a  careful  analytic  study  of  the  theory  of  resistive  interchange  modes  show  that 
there  are  several  very  important  stabilizing  effects  which  are  not  included  in  the  simple  slab  model  with 
gravity  simulating  curvature  and  have  been  ignored  in  previous  treatments  in  cylindrical  geometry. 

Most  resistive  interchange  modes  are  localized  to  within  the  singular  layer.  For  these  modes,  the 
effect  of  compressibility  and  shear  turn  out  to  be  strong  stabilizing  effects.  If  the  plasma  compression 
term,  proportional  to  V  p  (an  upper  case  V  is  used  to  avoid  confusion  with  the  growth  rate  y)  is  not 
negligible,  the  plasma  can  support  a  pressure  gradient  and  be  interchange  stable  even  in  the  absence  of 
shear.  This  was  first  pointed  out  by  Bernstein  et  al.4  for  a  plasma  with  Bz  -  0.  Also,  for  Bt  *  0  there 
is  some  residual  shear  stabilization  which  can  further  stabilize  interchange  modes  in  a  resistive  plasma. 
Thus,  most  of  the  resistive  interchange  modes  can  be  stabilized  by  this  combination  of  compressibility 
and  shear.  The  marginal  stability  condition  gives  rise  to  a  marginally  stable  relative  pressure  profile  in  a 
reversed  field  pinch. 
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However,  if  these  inodes  are  stable,  there  is  one  additional  unstable  mode  which  couples  to  the 
plasma  outside  of  the  singular  layer.  This  instability  exists  for  any  negative  pressure  gradients.  Above 
a  critical  pressure  gradient,  the  growth  rate  scales  as  y  —  like  a  conventional  resistive  inter¬ 

change  mode.  Below  the  critical  pressure  gradient  however,  y  —  (-p')4,  so  that  the  growth  rate  goes 
to  zero  extremely  quickly  as  p'  -*  0.  Furthermore  y  scales  as  |A'|-4,  where  A'  is  proportional  to  the 
jump  in  the  radial  derivative  of  the  perturbed  radial  magnetic  held,  so  that  short  wavelength  modes  are 
much  less  dangerous  than  the  longest  wavelength  (m  —  0  and  m  —  1)  modes.  For  all  practical  pur¬ 
poses,  this  critical  pressure  gradient  then  defines  a  marginal  stability  pressure  for  this  last  resistive  inter¬ 
change  mode. 

Section  II  reviews  the  equations  and  scalings  for  these  modes,  as  originally  derived  by  Coppi, 
Greene  &  Johnson  (CGJ).5  Section  III  derives  the  properties  of  the  localized  modes  and  calculates 
relative  pressure  profiles  for  a  reversed  field  pinch,  using  the  Bessel  function  equilibrium  model.  We 
conclude  that  the  marginally  stable  pressure  profile  can  have  a  central  pressure  up  to  400  times  the 
pressure  at  the  walls. 

In  Section  IV  we  discuss  the  properties  of  the  mode  that  couples  to  the  fluid  outside  the  singular 
layer.  In  several  limits  (high  shear,  low  shear),  assuming  parameters  so  that  the  previously  discussed 
localized  modes  are  stable,  the  dispersion  relation  takes  on  a  very  simple  form,  showing  the  y  —  (-p')4 
scaling  as  p'  -» 0.  We  also  show  that  the  dispersion  relation  can  be  obtained  very  simply  in  the  low 
shear,  low  pressure  gradient  limit.  We  evaluate  A'  for  modes  with  various  toroidal  and  poloidal  mode 
numbers  n.m,  using  the  Bessel  function  model  again  and  show  that,  if  the  pitch  parameter  is  large 
enough,  the  m  -  0  modes  are  the  most  dangerous.  This  result  was  observed  in  simulations  by  Schnack 
et  al.3  We  also  evaluate  the  critical  0  for  reversed  field  pinches.  The  critical  0  turns  out  to  be  about  7 
percent  for  a  magnetic  Reynolds  number  Rm  -  103  (T,  ~  10  eV)  but  lest*  than  one  percent  for  kilovolt 
temperatures  (J?„  >  10s).  The  former  result  agrees  with  the  numerical  results  of  Schnack  et  al.3  The 
critical  beta  decreases  with  conductivity  because  the  stabilizing  coupling  to  the  outer  (tearing  stable) 
region  becomes  weaker  with  increasing  conductivity. 
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II.  Review  of  the  Equations  A  Scalings 


In  this  section  we  briefly  review  the  Equations  and  scaling  as  derived  in  CGJ.  Outside  a  small  dis¬ 
tance  Lr  away  from  the  rational  surface,  the  ideal  MHD  equations  hold,  but  within  this  distance  I,, 
finite  resistivity  is  important.  As  we  will  see,  there  are  some  modes  which  are  localized  within  Lr  and 
which  do  not  couple  at  all  to  the  outer  region.  For  other  modes,  there  is  strong  coupling  to  the  outer 
region.  We  will  assume  that  Lr!a  —  e  where  a  is  the  radius  of  the  mode  rational  surface  and  a  is  also 
assumed  to  be  roughly  equal  to  the  scale  of  variation  of  the  eigenfunction  in  the  other  direction  perpen¬ 
dicular  to  B.  Since  the  plasma  is  assumed  to  be  ideal  MHD  stable  to  localized  modes  (Suydam  stable), 
the  growth  rate  also  scales  as  « .  In  this,  the  so-called  glow  interchange  ordering  of  CGJ,  €  also  scales  as 


In  cylindrical  geometry,  the  eigenfunction  has  the  form  f  (r)  expOyf  +  imO  +  i/cz)  where  in  a 
large  aspect  ratio  torus  k  —  -n/R.  A  superscript  —  denotes  a  perturbed  quantity.  The  mode  rational 
surface  is  determined  by 


m  *  nq 

where  q  -  rBjRB9.  The  operator  B  •  V  operating  on  a  scalar  is  given  by 


B  •  V  —  ■ 


in  B#  q* 


(r  -  a)  4-  0(r  -  a)2 


and  within  the  singular  region  is  of  order  e .  (We  use  this  form  rather  than  that  of  CGJ  since  it  applies 
also  to  m  -  0  modes,  which  have  their  rational  surface  where  Bz  -  0.) 

Since  the  resistive  interchange  mode  is  characterized  by  nonzero  resistivity  allowing  the  fluid  to 
slip  through  the  sheared  field  near  the  rational  surface,  the  mode  is  mostly  fluid  flow  with  only  a  small 
magnetic  perturbation.  Thus  B  —  c  |  where  |  is  the  linearized  fluid  displacement.  Since  the  growth 
rate  is  assumed  to  be  slow  compared  to  the  magnetosonic  speed,  we  have 


B  ■  B  +  p-  0. 


■  v-* 
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Since  p  -  -p'i ,  -  Vp  (V  •  { ),  where  T  is  the  specific  heat  ratio,  both  and  (V  •  ( )  must  be  of  order 
c .  Thus  to  zero  order  the  flow  is  divergence  free.  However  the  small  compression  does  play  an  impor¬ 
tant  role  in  the  pressure  perturbation.  Since  B  is  also  divergence  free,  B,  —  « 


If  we  define  a  local  coordinate  system  Jr,  ji  and  X  “X  x  ii*  where  j,  +ii  are  unit  vectors  in  the 
directions  of  r  and  £,  then 

£-lri-+flil  (4) 

M  ~  B,J,  +  Bw  Jt\  +  (5) 

where  |Hti  —  BL ,n  ~  «  and  B,  —  e 2.  The  equations  for  £,  B,  and  V  ■(  to  lowest  order  are 

given  by  CGJ  as 


+  /  R8‘  ~  -  0 
n  B  dx 

(6) 

RBe  d  Br 

*  +  i  nB  dx  -° 

(7) 

bb„-  p'ir  +  rp  (v  |) 

(8) 

,  „  -  .  i  n  Bgq'  Bx  - 

Py2B( -  p  B,  Bn 

a 

(9) 

0^i"r2nlBR  inB»<i'x  a- 
Py  *  '  R*a  8,1  a  B  ' 

(10) 

S  V  Bn  _  inB„q'x  . 

B'  yB'  a 

(11) 

_  2L  B"n  -  -  -■  B-q  -  -  B  (V  f) 

(12) 

~  (£  ±  B}Y  ; 

B 

where  x  “  r  -  a.  Equations  (6)-  (12)  are  CGJ  eqs.  (47)  •  (52)  with  several  slight  changes  in  notation: 
we  have  used  q  instead  of  the  rotational  transform  since  q  is  a  nonsingular  function  of  r  for  a  reversed 
field  pinch;  we  have  used  y  for  growth  rate  and  r  for  specific  heat  ratio;  finally  all  components  of  B 
and  |  have  dimension  magnetic  field  and  length  respectively. 

Equations  (6)  and  (7)  set  the  €°  part  of  V  •  |  and  V  ■  B  equal  to  zero.  Equation  (8)  sets  the 
perturbed  magnetic  plus  thermal  pressure  equal  to  zero.  Equation  (9)  is  the  dot  product  of  £  with  the 
momentum  equation.  Eq.  (10)  results  from  operating  on  the  momentum  equation  with  the  operator 
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V  *  ( B/B2x )  Coppi,  Greene  and  Johnson5  call  this  the  annihilated  momentum  equation  since  the 
operator  annihilated  information  about  propagation  along  the  lines  of  force  (given  by  eq.  (9))  and  also 
annihilates  information  about  V  (B  B  +  p)  (given  by  eq.  (3)).  Equations  (ll)  and  (12)  are  the 
radial  and  ji  components  of  Ohm's  law. 


The  first  two  terms  of  (12)  combine  to  give  -B  (Vx  •  fA),  the  compression  of  the  field  due  to  the 
compression  of  the  fluid  perpendicular  to  B.  The  last  term  is  a  combination  of  two  effects.  First  there 
is  the  convection  of  B  into  regions  of  different  B.  This  gives  a  contribution  to  B\\  of 

|ll|.  (13, 

Secondly,  consider  a  flux  tube  at  radius  in  which  is  incompressibly  displaced  a  distance  |r.  If  the  tube 

A  - 

originally  has  area  A ,  the  change  in  area  of  the  displaced  tube  is  8  A  —  — -  f  r.  However  since  the 


flux  is  conserved,  then 


a  ft  -  -  4  •  A  “  “  <r 

A  r 


_  £r  j|_  .  B* 

B  Br  P+  2 


The  last  relation  follows  from  pressure  balance.  Combining  Eqs.  (13)  and  (15)  gives  the  last  term  in 
Eq.  (12).  The  term  on  the  right  of  (11)  is  the  analog  of  the  first  term  on  the  right  in  (12)  or  ,  alterna¬ 
tively,  the  convection  of  the  helical  flux  4/  =  t  r  Br. 

Equations  (8)  and  (9)  give  V  •  |  and  In  in  terms  of  the  other  variables.  In  the  notation  of  CGJ, 
the  equations  are 

QW  +  X  S)  (16) 

+  «7, 


Y ^l1  +  +  nf) Y  +  ° -  fj) a  +%*■* 


where  Q  -  y/Qr,  H  -  f  „  ¥  =  -  ^~,L 


B„X-x/L„  Y-  jMh  «n.«nd 
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D _ 

Bt2R2q2 


(19) 


fi  -  2 p(a)/B2 
S-  */R2q‘2 


(20) 

(21) 


Lr- 


n2  Bj  q  ‘ 


1/6 


Qr~ 


tj  n2  B»  q2 


1 1/3 


(22) 

(23) 


P  a 


We  have  eliminated  m,  Bz  and  i  —  2 n/q  because  in  a  reversed  field  pinch  the  mode  rational  surface  for 
m  -  0  modes  has  Bz  —  0. 


As  long  as  there  is  no  coupling  to  the  outer  region,  the  solutions  to  Eqs.  (16)  -  (18)  have  either 
odd  or  even  symmetry,  with  the  symmetry  of  ♦  opposite  to  the  symmetry  of  S  and  Y. 


As  shown  in  CGJ,  near  marginal  stability  these  modes  are  localized  within  a  more  narrow  layer, 
with  X  —  Qv*,  and  therefore  modes  with  odd  S,  even  ♦  and  Y  in  fact  must  be  treated  by  tearing  ord¬ 
ering.  The  most  important  effects  of  this  ordering,  where  y  —  t>3/5,  r  -  rs  —  r)v$  and  /3  —  is  that 
the  first  term  on  the  right  in  (16)  can  be  treated  to  lowest  order  as  a  constant.  (This  is  the  constant-^ 
approximation  of  Ref.  6.)  Also,  the  first  term  in  the  Y  coefficient  in  (18)  and  the  second  term  in  the  E 
coefficient  are  negligible. 


III.  Even  E  Modes  and  Stable  Pressure  Profiles 

In  this  section  we  review  the  even  E  modes.5-6  For  modes  of  this  symmetry  it  is  consistent  with 
the  constant  ¥  approximation  to  neglect  the  first  term  of  (11)  i.e.,  the  first  term  on  the  right  in  (16). 
Therefore,  Eqs.  (17)  and  (18)  decouple  from  Eq.  (16)  for  ¥.  With  tearing  ordering  (i.e.,  B  ~  tj275), 
Eqs.  (17)  and  (18)  become 
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Following  CGJ  we  expand  in  Hermite  functions  3  -  £3,  H,(z)  exp(-zJ/2),  Y  -  £  Y,  H,(z) 

i  i 

exp  (-  z2/ 2),  where  z  -  X! Q'/4. 


Making  use  of  the  fact  that 


one  finds  for  even  /  two  simultaneous  equations  for  3/andY,, 


H,(z)  exp  -  -y 


-  -  (2/  +  1)  Ht(z)  exp  (-  z2/ 2) 


-(21  +  1)  5,  +  Q-vtY,-  0 


!|s  -  3/  + 

-(21  +  1)  - 

r/9  J 

r/s 

Y,-  0 


(26) 


which  give  the  dispersion  relation 

It 

v  2 

The  most  unstab/e  mode  has  /  —  0;  therefore  stabthy  to  all  even  S  modes  is  guaranteed  if 

I?  <  5  +  1  (28a) 


20 

It 


-  S 


2/  +  l 


-  (2/  +  1) 


(27) 


! 


or  in  physical  units 


_ *£>_  <  l£ 

/?2fl,y2  a2 


(28b) 


The  even  (,  modes  have  I,  and  B\\  entirely  localized  to  a  region  X  <  Ql/*  inside  the  singular  layer  and 
do  not  couple  to  the  plasma  in  the  outer  region.  To  check  the  validity  of  the  constant  ¥  approximation 
we  need  only  show  a  postiori  that  Qi>  <<<!>".  Using  —  'V/QU1,  we  find  the  condition  for  validity 
of  the  constant  approximation  is 


Qi/2  —  d  —  0  «  1 


(29) 


Notice  that  unlike  a  gravity  driven  resistive  interchange  in  slab  geometry6,  the  even  3  modes  in  a 
cylinder  are  stable  for  a  pressure  gradient  below  a  critical  value  given  by  Eq.  (28b). 
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This  critical  pressure  gradient  can  be  related  to  previous  work4  on  the  ideal  MHD  stability  of  a 
compressible  diffuse  pinch  with  Bz  -  0.  According  to  Bernstein  et  al.\  a  diffuse  z  pinch,  which  has 
unfavorable  curvature  everywhere,  is  stable  to  a  pure  interchange  mode  only  if 

Ir  +  {~p  >  0  (30) 

where  V  is  proportional  to  the  flux  tube  volume  r/B9.  Thus,  increasing  T  p  (i.e.,  compressibility)  has  a 
stabilizing  effect.  It  is  not  difficult  to  show  that  (30)  is  proportional  to  the  factor  S  -  D( 2/H3  +  1)  of 
Eq.  (18)  if  Bz  -  0  is  assumed.  The  third  term  in  this  factor  is  necessary  because  0  must  be  of  order 
unity  if  B2  *  0. 

Increasing  Bz  has  both  stabilizing  and  destabilizing  effects  on  even  5  modes.  The  destablizing 
effect,  given  by  the  factor  T  p/  B2  in  (28b),  is  due  to  the  fact  that  large  B2  implies  incompressibility  (but 
without  C 2  -  T  p/p  —  <»)  so  that  no  work  is  done  in  compressing  the  plasma.  The  stabilizing  effect  of 
increasing  Bz  is  that  this  may  increase  the  shear  [represented  by  the  second  term  on  the  right  of  Eq. 
(28a)  or  Eq.  (28b)]. 

To  get  an  idea  of  the  stable  pressure  profiles  in  a  reversed  field  pinch,  we  have  computed  the  mar¬ 
ginal  pressure  profile  from  Eq.  (28b)  for  a  Bessel  function  model 

B0Jobr)  (31) 

Bq  *  Bq  r). 

This  is  the  magnetic  field  produced  by  the  plasma  if  it  relaxes  to  minimum  energy  state  while  maintain¬ 
ing  constant  helicity7,  giving  V  x  B  -  p  B.  Although  the  field  given  by  Eq.  (31)  is  force  free,  it 
should  still  be  a  good  approximation  to  the  field  structure  in  a  low  0  plasma.  If  the  quantity  rH  is  the 
wall  radius,  the  range  of  m  of  most  interest  to  a  reversed  field  pinch  is 

2,4  <  p.  rw  <  3.1, 

the  former  value  being  the  value  required  for  field  reversal,  the  latter  required  for  tearing  mode  stabil¬ 
ity.8 


8 


.  n  Vfc^Sin-— "r 
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Since  Eq.  (28b)  gives  pip  it  gives  only  relative,  not  absolute,  pressure  profiles.  Thus,  if  the 
boundary  condition  is  p  (/*„)  *  0,  there  must  be  an  unstable  region  somewhere  in  the  plasma.  In  Fig. 
(1)  are  shown  plots  of  relative  pressure  profiles  for  five  different  values  of  m  Clearly  as^  increases, 
the  pressure  at  the  center  can  become  quite  large  compared  to  the  pressure  at  the  walls.  For  instance  at 
M  *  31,  just  below  the  threshold  for  tearing  mode  instability,  the  pressure  at  the  center  is  more  than 
400  times  the  pressure  at  the  wall.  It  is  of  interest  to  note  that  the  toroidal  factor  H  of  Ref.  9  [eq. 
(13)1  is  zero  for  toroidal  equilibria  which  satisfy  V  x  B  -  m  B  with  p,  -  j  *  B/ B2  constant. 

IV.  Odd  E  Modes  — The  Maximum  Stable  Pressure 

In  the  previous  section  we  have  seen  that  if  the  pressure  gradient  is  smaller  than  that  given  by  Eq. 
(28a  or  b),  then  the  even  E  modes  are  stable.  One  might  also  think  odd  E  modes  stable  also  since 
they  have  larger  /.  However,  since  S  has  odd  symmetry,  V  has  even  symmetry,  so  it  can  couple  to  the 
outer  MHD  stable  region  through  the  (constant)  value  of  This  allows  for  the  presence  of  one 
unstable  mode  even  if  (b)  is  violated. 

The  inner  region  equations  for  odd  S  modes  are 

S'  -  z2 E  +  Q-}/2 Y  -  0'1/4z^o  (32) 

Y"  -  z2 Y  -  QV2  (S  ~  2£>/r/3)  E 
-  2Qyi\/Tp  -  DQ~u*z  ¥<,.  (33) 

Henceforth,  we  assume  that  the  even  E  modes  are  stable,  i.e.,  S-  2D/T0  >  -  1.  Earlier  refer¬ 
ences5  6  have  shown  that  this  equation  can  be  analyzed  by  expanding  in  Hermite  functions  and  express¬ 
ing  the  resulting  series  as  a  hypergeometric  series.  In  the  incompressible  limit  2D/V  0  —  0,  we  find,  as 
in  Refs.  5,  9, 

L'i  -  TriraW  w> 

where  A  -  (d'P/dr  (<j+)  -  (a—)]/'V  (a).  In  the  limit  S  »  1  (in  which  2D/T 0  need  not  be 

small,  and  which  corresponds  to  the  low  shear  limit)  we  find 
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^  (Qs,*cr 1/4  —  D/QV*(t ^4), 

where  a  -  S  -  2D/T0. 


(35) 


In  the  Appendix  we  show  that  it  is  possible  to  find  a  simple  closed-form  solution  to  (32),  (33)  in 
this  latter  limit.  In  fact,  we  find 


H  (z) 


VqQ~1/4z(D/Q>/2  +  z2) 
z4  +  a 


(36) 


We  are  also  able  to  obtain  (35)  exactly. 


It  has  been  shown  by  using  Nyquist  analysis  that  equations  of  the  form  (34)  or  (35)  always  have 
one  root  in  the  right  half  Q  plane.9  A  simple  graphical  analysis  shows  that  this  root  is  pure  real.  If 
A'  >  0,  D  «  1,  the  A'  balances  the  05/4  term  on  the  right  hand  side  of  eq.  (34)  to  give  the  conven¬ 
tional  tearing  mode  dispersion  relation.  On  the  other  hand,  if  A'  -  0,  one  easily  calculates  from  (34) 
Q  -  (tr  D/ 4)2/3,  the  conventional  result  for  a  resistive  interchange  mode.  Using  (35)  we  obtain  a  simi¬ 
lar  result  Q  -  Dv*/(rvy.  The  most  interesting  limit,  however,  is  for  A'  negative  and  either  |A'|  »  I 
or  D  «  1.  In  this  case  (34)  gives 

Q  7.1  (D/ LAV.  (37) 

This  growth  rate  goes  to  zero  very  rapidly  as  D  — *  0.  The  result  complements  that  of  Glasser,  Greene 
A  Johnson9,  who  show  that  coupling  to  favorable  average  curvature  can  stabilize  a  tearing  mode.  Here 
we  show  the  converse,  that  coupling  to  large,  negative  A'  can  substantially  reduce  the  growth  rate  of  a 
resistive  interchange  mode.  The  growth  rate  goes  to  zero  so  rapidly  with  D  that  we  can  consider  the 
transition  point  from  Dv*  to  D4  behavior,  namely  Dc  -  .52  (—  Lr  A')^5,  to  be  an  effective  critical  beta 
for  the  mode.  Using  (35)  the  corresponding  results  are  Q  =*  (24/cr)  (D/ LAV  and 
Dc  -  .3&r 1/5  (—  Lt  A')6^5.  Note  that  (37)  gives  y  —  p'**!  whereas  the  <r  —  <»  result  gives 
y  ~~  P  4I v  Q  4-  In  both  cases  y  —  tj  3/5 . 
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in  physical  variables  the  critical  pressure  gradient  obtained  from  Dc  is 

-  dp! dr  -  (Bi  R2q2/ 4)  (-  A'L,)675  (38a) 

2  2  ),/5 

-  <-A')*5  -^L«  .  (38b) 

n  q  B§  J 

Note  that,  for  plasma  near  the  critical  0  for  low  m.n  modes,  modes  with  larger  m,n  have  a  much  smaller 
growth  rate  since  Q  scales  as  (A')“4  and  since  A'  —  —  n  for  n  —• ►  «>.  in  particular,  since  Q,  —  n273,  we 
have  y  —  *“1Q73  as  n  — - ►  Another  way  of  expressing  this  is  to  note  that  the  critical  pressure  gradient 

scales  as  /i475  for  .large  n.  Note  that  the  critical  pressure  gradient  scales  as  tj275.  This  is  as  expected, 
since  it  is  tearing  ordering  with/3  —  t?Z  5  that  makes  the  three  terms  of  (34)  or  of  (35)  comparable. 

For  the  force  free  Bessel  function  model  B9^BQJxfar),  Bz «  B0J0  ip  r),  we  can  easily 
compute  A'  analytically.  This  is  possible  because  the  outer  region  equations  for  this  equilibrium  are 
merely  V  x  B  -  p.  2?,  with  the  same  value  of  m  as  the  equilibrium.  (Recall  from  Sec.  Ill  that  toroidal 
generalizations  of  such  equilibria  have  the  toroidal  term  H  of  Ref.  9  equal  to  zero.)  In  this  mode!  field 
reversal  occurs  if  p.  rw  >  2.4  (rw  is  the  plasma  minor  radius)  and  an  m  *  1  tearing  mode  is  unstable  if 
p  rw  >  3.1.  In  Fig.  2a  we  show  the  maximum  value  of  A'r*  (  over  0  <  m  <  5,  -25  ^  n  ^  25)  as  a 
function  of  r  for  m  rw  «  2.81.  Fig.  2b  shows  the  maximum  of  A'r*.  over  n  for  various  m,  as  a  function 
of  p  rw.  These  results  show  that  for  2.4  <  p  rH  <  2.83,  the  m  ~  1  mode  with  n  >  0  (i.e.,  with  mode 
rational  surface  inside  the  field  reversal  point  )  has  the  largest  A'.  The  m  -  I  mode  with  largest  A'  has 
n  rjR  =  2,  which  has  the  effect  of  putting  the  mode  rational  surface  as  far  from  the  walls  and  from 
the  origin  r  —  0  as  possible.  For  2.83  <  p  rw  <  3.1,  however,  the  m  -  0  mode  (which  has  its  mode 
rational  surface  at  the  field  null  p  a  —  2.4  regardless  of  n )  with  n  *  ±  1  has  the  largest  value  of  A'. 
We  also  observe  that  an  m  —  1  tearing  mode  with  n  <  0  (i.e.,  with  mode  rational  surface  outside  the 
field  null)  becomes  unstable  (i.e..  A'  >  0)  for**  rw  >  3.1,  in  agreement  with  previous  studies.1  8 

Putting  A'1*  -4,  we  find  the  critical  Dfor  m  ■»  0  modes  to  be  approximately  Dc  -  3  fl'  275,  where 
Rm  is  the  magnetic  Reynolds  number.  For  Rm  -  10\  as  in  Ref.  3,  we  find  Df  —  20  percent.  For 
Rm  -  105,  we  find  Dc  =  3  percent.  For  the  marginal  equilibria  of  Sec.  Ill,  we  have  0  —  D/3,  giving 
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critical  0  equal  to  7  and  one  percent,  respectively.  The  former  result  is  in  agreement  with  the  simula¬ 
tions  of  Schnack  et  al.3 

In  summary,  we  find  that  matching  to  a  tearing  stable  exterior  region  (i.e.,  with  A'  <  0)  decreases 
growth  rates  for  the  resistive  interchange  sufficiently  that  below  a  critical  value  of  the  Suydam  parame¬ 
ter  A  the  growth  rate  scales  as  D 4  and  is  therefore  quite  small.  We  also  find  that  in  this  regime,  the 
growth  rate  y  scales  as  |A'|"4,  so  that  for  very  short  wavelength  modes  (n  — • >  »),  y  scales  as  n"10/3. 
From  a  calculation  of  A'  from  the  outer  region,  we  find  that  the  most  unstable  mode  is  either  an  m  -  1 
mode  with  n  rj R  =  2  or  an  m  —  0  mode  with  n  —  ±  1.  The  above  results,  notably  the  abrupt 
dropoff  in  growth  rate  for  small  A  the  decrease  in  growth  rate  for  large  n  and  the  predominance  of  the 
m  -  0,  m  -  1  modes,  are  in  agreement  with  the  numerical  simulations  of  Schnack  et  al.3 
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From  this  we  obtain 

Lr  A  -  J  dz  [♦o  +  Ol/4zH  (r)|<fc  (A7) 

or 

LA  “  72  1/4  <?V4  - 1>/»  ,/4o1/4| .  ( A8) 

which  is  exactly  (35).  For  S  >  >  1  and  2D/?P-\,  Qin  is  bounded  by  D/Sin,  so  the  third  and 
fourth  terms  of  (33)  are  in  relation  S:  Qin/ D  l  using  (32)]  or  S  :S~U1.  This  proves  that  the  approxi¬ 
mation  discussed  in  the  first  paragraph  of  this  appendix  is  valid. 
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Fig.  I  —  Marginally  stable  relative  pressure  profiles  (even  5  modes)  Tor  the  Bessel  function  equilibria 
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Fig.  2b  —  Maximum  value  of  A’  rw  over  n  for  m  “  0,  m  "*  1  modes  with  €  =  rHJ  R  Notice  that  tri  —  0 

modes  have  the  larger  A'  Tor  ft  rw  >  2.83 
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